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ABSTRACT 

As observational evidence increasingly consolidates the case for a cosmological constant A 
being the source of the Universe's accelerated expansion, the question whether, and if so, 
how well, future experiments could detect deviations from this standard scenario is raised 
with urgency. Assuming a dark energy component different from a cosmological constant, the 
observable effects in general include gravitational clustering described by the fluid's (rest- 
frame) speed of sound c,. We employ 3d weak cosmic shear, a proposed method to take 
advantage of the full three-dimensional information inherent to the cosmic shear field, to 
explore the capability of future surveys to detect dark energy clustering and the signature of 
an enhanced amplitude of the matter power spectrum on large scales. For this purpose, we 
present adequate numerical methods facilitating 3d weak cosmic shear calculations. We find 
that the possible constraints heavily depend on the dark energy equation of state w. If w is not 
very close to -1, constraining the squared sound speed c], within an order of magnitude seems 
possible with a combination of Euclid and Planck data. 
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1 INTRODUCTION 

Explaining the accelerated expansion of the Universe is one of the 
key tasks of cosmology today. If Einstein's general relativity re- 
mains unaltered on cosmological scales, the observed accelerated 
expansion — if no local effect — is due to a cosmological constant 
or an unknown cosmological fluid with negative pressure, the dark 
energy. To this day, all major observations are consistent with a 
cosmological constant jKomatsuetai]|201ll : iBartelmanijboiOah . 
Its unexpectedly tiny value — the cosmological constant problem — 
and the fact that its energy density is comparable t o that of mat - 
ter just today — the coincidence problem — (see, e. g.. lCarrolll200lh 
motivate the search for alternative models of dynamically evolving 
dark energy. 

Due to the lack of observational evidence for inhomogeneities 
in the dark energy, most studies have only investigated the conse- 
quences of a perfectly homogeneous dark energy component. Such 
a fluid is completely determined by its energy density and its equa- 
tion of state w. Its direct efi^ect is restricted to the expansion history, 
which indirectly causes a scale-independent modification of the 
growth rate of matter perturbations. In general, however, a cosmo- 
logical fluid can also leave signatures, possibly scale-dependent, 
by virtue of its perturbations. This could be, in principle, a means 
to discriminate between different dark energy models. 

Once we include linear perturbations, another characteristic 
quantity enters the scene, the sound speed c,. This quantity de- 
fines a sound horizon such that scales outside and inside this hori- 
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zon can undergo different evolutions. In general, both w and c, are 
necessary to describe the observable effects of a fluid. In order to 
explore the na ture of the dark energy, cosmolog y has to constrain 
both w and c, jHull2002bl : lErickson et al. 2002; DeDeo et al.ll2003l : 



, p 

IHu & Scran"toi]|2004l) 

The remarkable progress of observational cosmology in de- 
termining the fundamental parameters describing our Universe has 
not yet led to significant constraints on the dark energy sound speed 
c, ISean&Dore 2004; de Putter et al. 2010; Li & Xia 2010). In 
this work, we study whether next generation precision observations 
of the cosmic microwave back ground together with the proposed 
method 3d weak cosmic shear iHeavensll2003h have the potential 
of providing significant progress in this respect. 

3d weak cosmic shear is a metho d to gain precision informa- 
tion a b out the growth of pert urbations jHeaven sll2003l ;l Castro et al.l 
|2005| ; iKitching et all I20T J). Contrary to ordinary galaxy sur- 
veys, it has — like weak lensing in general — the advantage of be- 
ing independent of galaxy bias models. Only well-understood 
general relativity is needed from the theoretical side. This is 
one reason why weak gra vitational lensing, since its beginnings 
jVan Waerb eke et al."2000l). has advanced along with the CMB to 
one o f the co smologi cal p robes with the largest potential ( HutereJ 
l20ld ; I BartelmannI l2010bl) . Weak lensing methods have in fact 
proved to be powerful tools to constrain dark energy, i. e. mainly 
its equation of state parameter w l^ uterer 2002; Jain & Tavlog 
'2003'; Heavens 2003; Bernstein & Jain 2004; Takada & Jain 200l 
Hannestad et al. 2006; Heavens et al. 2006; Amendola et al. 200g; 
IlloUenstein et al.i2009l ; iKilbinger et alj|2009i ; lHutereill20ia) . 

Most weak lensing studies consider the case of tomographic 
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measurements where the sample of lensed galaxies is split up into 
redshi ft bins on which t he standard weak lensing methods are ap- 
plied i IhuII1999L 120023) . The advantage of tomography is an en- 
hanced sensitivity due to reduced averaging along a line of sight 
compared to unbinned cosmic shear spectra, but the shape of the 
dark matter power spectrum is not measured independently from 
growth factors and geometry (recent studi es about tomography 
and the relation to 3d weak lensing include iKitching et alj|2011l : 
ISchafer & Heisenberg|l201 ll) . 

The 3d version of weak lensing is a complement to standard 2d 
weak lensing with the aim of retaining the full three-dimensional 
information contained in the cosmological shear field. The starting 
point is to not only make use of the angular positions of lensed 
galaxies on the sky, but to also include their redshifts as a distance 
measure such that each individual galaxy provides a measure of the 
tidal shear. 

Let us briefly explain our motivation to look specifically into 
3d weak lensing as opposed to tomographic methods. Weak lensing 
spectra provide an integral measure of the dark matter power spec- 
trum, weighted with the lensing efficiency function. The enhance- 
ment of the matter power spectrum due to the clustering of dark 
energy is restricted to large scales and would thus influence a weak 
lensing convergence spectrum only little. A 3d method, however, 
provides a direct measurement of the amplitude of the dark matter 
spectrum and would be better suited to distinguish enhanced spec- 
tra from unenhanced spectra and therefore to provide constraints 
on the properties of dark energy and its clustering. This would ef- 
fectively break the degeneracy between the power spectrum shape 
and the lensing efficiency, consisting of the growth function and 
geometrical factors, such that the signature of dark energy induced 
clustering should be easier to observe. 

The organisation of this paper is as follows. We first describe 
clustering dark energy in general, make contact with prominent 
dark energy models, and introduce a parameterisation in Sec.|2] We 
then explain the 3d weak cosmic shear method in Sec. |3] A brief 
description of the Fisher matrix method for forecasting parameter 
constraints is given in Sec.|4] We present adequate and efficient nu- 
merical tools in Sec. [5] Our results are shown in Sec. [6] and we 
conclude in Sec.|7] 



2 CLUSTERING DARK ENERGY 
2.1 The sound speed 

The dynamics of the background and the evolution of scalar lin- 
ear perturbations of a cosmological fluid are fully determined by 
its equation of state w = pip and its (squared) sound speed 
c] = 6p/Sp. If we describe dark energy as a cosmological fluid, 
coupled to other fluids only by virtue of the gravitational interac- 
tion, the natural parameters are w and cj. 

The sound speed defines a characteristic scale A oc |c.,|, be- 
low which the fluid resists gravitational collapse. In turn, this means 
that the effects of gravitational clustering are only observable if the 
scale A lies within the Hubble horizon, A < H^' , where H = a' /a 
is the conformal Hubble parameter and a prime denotes a derivative 
with respect to conformal time t. 

In general, the speed of sound is defined by the quotient of 
the pressure and density perturbations, = 6p/dp. Both, 6p and 
dp, however, are gauge-dependent quantities, whence we shall only 
consider the gauge-invariant rest-frame speed of sound defined in 
a frame where the velocity perturbation of the fluid vanishes, v = 0. 



We can illustrate the role of the sound speed with the help 
of the evolution equations of linear perturbations, which are ob- 
tained from the general energy-momentum conservation equations 
T''^-y = 0. These equations are valid if there is no coupling, i. e. no 
energy-momentum exchange, between the fluid and other compo- 
nents such as matter. As usual, we split into background quantities 
and linear perturbations, T'"' = f*"' -I- 6T''^, and we define T^o = 
-p(l + 6), T°i = (p + p)vi, T'j = {p + 6p)5) + S'y. We further define 
a gauge-invariant density perturbation A = 5 + 3(1 + w)y(l' - S) 
in Fourier space, where B is a metric perturbation defined as in 
iKodama & Sasakil (Il984h . Choosing the fluid's rest frame v = Q 
and a gauge where fi = 0, we simply get A = 5. Describing the 
evolution of perturbations, for a single fluid, in terms of the vari- 
able pa^6 , we find the following se cond-order differential equation 
(cf., e. g.. lKodama & Sasakilll984) , 

{pa'S)" + {i + 'i^'H(pa^S)' + {lc-c] - ^(l+wYH^ipa^S) = 0,(1) 

neglecting anisotropic shear, E'^ = 0. A critical scale ^crit = 1 1'^au 
is given by the vanishing of the source term oc pa^S driving gravi- 
tational collapse, i. e.. 

The perturbation variable pa^6 can only grow on subhorizon scales 
for /icrii < "H"', which translates into the approximate relation c] < 
i + w. Especially for an equation of state w close to - 1, as preferred 
by current observations (Komatsu et al. 2011), this only occurs for 
very small sound speeds <c 1. These effects are restricted to 
large scales A > A^^n. 

In a complete description of the perturbation evolution, we 
have to cope with the multi-component fluid of (at least) matter 
and dark energy. Nonetheless, we can still motivate a corresponding 
heuristic definition of an effective scale characterising dark energy 
clustering, see Sec. 12.31 We will then also show quantitatively how 
a clustering dark energy component (with constant w and c^) affects 
the large-scale matter power spectrum P(k). 

The (rest-frame) sound speed considered here may not be 
confused with the adiabatic sound speed c^, which is only equal 
to the quotient Sp/6p for adiabatic perturbations, i. e. when the en- 
tropy perturbation is zero. In general, it is given by cl = pip. The 
difference between the two quantities defines a gauge-invariant en- 
tropy perturbation {c] - cl)6/w. For a fluid with constant equation 
of state w, the adiabatic sound speed simply reduces to cl = w. For 
a brief introduction to dark energy clustering, see lOordon & Hul 



2.2 Relation to common dark energy models 

2.2.1 Quintessence 

The most pro minent example of dynamical dark energ y is standard 
quintessence l lWettericri[T98^ : iRatra & Peeblesi f 198m), i. e. a cos- 
mological scalar field (p with standard kinetic term and a potential 
V(ip), defining a Lagrangian density £ = -^dfiifid^'if - V(ip). For 
suitable choices of the potential V(ifi), the dynamics of the back- 
ground field ip shows appealing tracker behaviours providing ro- 
bustness against initial conditions. 

The perturbation S(p of the quintessence field usually is of 
little importance on subhorizon scales, the reason being that the 
quintessence sound speed cj is unity. 
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This is easily seen by explicitly writing energy density and 
pressure perturbations of the scalar field, 



6p^ = tfiStp - V_^S<p. 



(3) 
(4) 



Since the velocity perturbation v is proportional to the field per- 
turbation Sif, the rest-frame speed of sound (for v = Stp = 0) is 



1. 



We conclude that the detection of a dark energy sound speed 
< 1 would not only challenge the ACDM model but standard 
quintessence models as well. 

A class of models with very diffe rent behaviour, however, is 
given by coupled quintessen ce models jWetterichl 19951 ; lAmendolal 
llOOdlAmendoia et al.ll2008[) . In these models, there is an energy- 
momentum exchange between the dark energy and other compo- 
nents such as dark matter or neutrinos. Dark energy can then no 
longer be described as an independent fluid, and the equations 
of Sec. 12.11 do not apply. In fact, subhorizon perturbations of the 
quintessence field can grow in these models. Although not consid- 
ered in this work, the case of energy-momentum exchange between 
dark energy and matter has been parametrised and studied in the 
light of weak lensing (Schiifer et al. 2 008 ; La Vacc a & C olpmbo 
I2OO8I ; ICaldera-CabraleralT2009l ; jPe Bemardis et al.ll201 ih . 



2.2.2 k— essence 

Looking at Eqs. ([3]l and (|4j, the reason for = 1 in standard 
quintessence is the identical dependence of 6p^ and Sp^ on Sip. For- 
mally, this could easily be changed by allowing the potential to also 
depend on i^, V = ip). If this dependence can be split into two 
sunmiands, we could reinterpret the (p dependence as a modifica- 
tion not of the potential but of the kinetic term. 

Non-standard kinetic terms are th e starting point for k- 
essence mo dels of dynamical dark energy jArmendariz-Picon et al.l 
I2000ll200lh . In these models, the Lagrangian X is a generic func- 
tion of the standard kinetic term X = -^d^ip&'ip. It is thus possible 
for the sound speed cP, to take any value, without violating causality 
teabichev et"al]|2008l) . 

The energy density and the pressure are given by the corre- 
sponding components of the energy-momentum tensor. They read 
p = IJLjiX - £. and p = £.. The equ ation of state w = pi p and the 
rest-frame sound speed = Sp/dp terickson et alj|2002l) are then 



£ 



£x 

£,x + 2£,xxX' 



(5) 
(6) 



Of course, both w and evolve in time and may take very different 
values at different epochs. The question whether the time evolution 
of c] could leav e characteristic observational imprints was studied 
by lAnsari & U nnikrishnan ( 2011). In the framework of a specific 
fc-essence type model, 3d cosmic shear has b een used to forecas t 
possible constraints on the model parameters dCamera et al.ll2O10h . 



2.3 Parametrised clustering dark energy 

A frequently employed parameterisation of dark energy, which we 
shall adopt here, is the wCDM model, sometimes called XCDM 
( lTumer&Whit3ll997h . Contrary to a cosmological constant A 
with equation of state Wa = -1, the model allows for an arbitrary 
dark energy equation of state w which is taken to be constant in 



time. The model is often extended to a linear evolution of w with 
respect to th e scale factor a (for an attempt to study w as a free 
function, cf. iHuterer & Turneill200ll) . The simplest generalisation 
for including possible clustering of dark energy is to further intro- 
duce a rest-frame sound speed c], also constant in time. In this 
paper, we completely parametrise the dark energy component by 
constant numbers w and Cj. 

Dynamical dark energy such as quintessence and fc-essence 
provides a large class of models that cannot be approximated by 
a simple parameterisation such as the wCDM model. In fact, the 
wCDM model (for constant w) does not resemble very closely any 
of the prominent dynamical models. Whenever new observational 
data are published, it is thus not sufficient to study constraints in the 
vvCDM model alone, but to also study the individual dark energy 
models. 

Nonetheless, the wCDM model is, in terms of its parameters, 
a somewhat minimal extension of the standard ACDM model, in- 
cluding the latter as a special case. Hence, it is a useful tool to 
forecast how strong the deviations from ACDM must be for future 
observations to detect them. 

We now turn to the description of linear perturbations in 
the presence of a clustering dark energy component parametrised 
as above. The linear growth of perturbations is described by a 
growth function g{k, a) that links the Newtonian gravitational po- 
tential <l>t(a) at scale factor a to the one today according to 
Ot(a) = g(k,a)(t>'j^/a. We shall now introduce the parameterisation 
for g(k, a) that we use for the study of dark energy perturbations in 
the wCDM model. 

Dark energy perturbations contribute to the gravitational po- 
tential just as matter perturbations via the Poisson equation, 

yt^O = -4nGa^ (p™A<'"> + Pde^^"'^'') = -4nGa- ep,„A<"'', (7) 

where we have used the gauge-invariant density perturbations A*""' 
and A*"^' and introduced the quantity Q = Q(k, a). It is defined via 



2=1 + 



Pde 



(8) 



p„,A<'") 

An important effect of dark energy perturbations is their influ- 
ence on the growth of matter perturbations expressed in terms 
of a modified growth index 7 defined by dlnA*"'Vdlna = Q.J„ 
jLinder & C ahn 2007j)^ As a function of Q, we may approxi- 
mate (cf. lSapone & Kunzl2009l ; ISaDone et al.l2010l ; rLinder & Cahi] 
l2007b 



r 



3(1 -w-A) 



, A ■ 



e-1 



(9) 



5 - 6w 1 - n„ 

We follow ISapone et al.l ( l2010h . parametrising Q(k,a) for the 
wCDM model with sound speed c\ as 

2 ecla 



a , i-o!;,(i+w)fl-^"' , 



3 n" ' 

'^m ' '0 



(10) 



Together, these equations provide a convenient way of obtaining 
the growth function 



g(fc,«) = ^^exp( I r:- «„,(«')''**■"■' I . (11) 

The explicit appearance of Q in this expression is due to our defi- 
nition of g describing the growth of the total gravitational potential 
rather than of the matter perturbations only. 

Of course, the growth function could easily be directly ob- 
tained by solving the linear perturbation equations numerically. For 
illustration, we show the linear matter power spectrum P(K) for the 
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Figure 1. The upper figure shows the matter power spectrum P(k) for w = 
-0.8 and varying dark energy speed of sound cj divided by the spectrum for 
= 1 . The matter power spectrum is computed from the gauge-invariant 
density contrast A''"'. The lower figure shows the scale dependence of 
for the same values of w and cl. 



wCDM model for different sound speeds and w = -0.8 in Fig.[T] 
Here, we have used the code camb (Code for Anisotropies in the 
Microwave Background, Lewis et alJflOOO '), which has built-in fa- 
cilities to work with the wCDM model with constant c^. We have 
assumed adiabatic initial conditions. 

Perturbations in the dark energy act as an extra source of 
the gravitational potential in the Poisson equation ^ enhancing 
the growth of matter perturbations on subhorizon scales. This en- 
hancement, however, is less than a percent effect for sound speeds 
> 0.1 and restricted to large scales. The power spectrum P(k) 
has to be known with very high precision in order to find significant 
constraints on c^. This becomes even more difficult for w closer to 
-1, cf. Eq. l2j. Note that the plot also shows superhorizon scales, 
where the results are gauge-dependent. For our analysis, we will 
use (subhorizon) scales A: between 10"^ and 10"' Mpc"'. 

With the parameterisation of Q(k, a), Eq. HOI , at hand, we 
can ask above which scale A^ff dark energy clustering could leave 
observable traces. Let us make the heuristic assumption that the 
effect of a clustering dark energy component would be observable 
once roughly Q(k, a = 1) > 1 + £, with e for example at the percent 
level. This is the case for scales 



1 / £ 



1/2 



1 



3 VI +w'Mq 



(12) 



For the exemplary values e 
critical scale 



'Hq- 



« 0.1 



vr 



1% and nl a 0.3, this defines a 



(13) 



with a similar behaviour as the scale given in Eq. A precision 
experiment might be able to detect dark energy clustering if the 
effective scale A^g lies within the Hubble horizon. In particular, the 
common choice w = -1 refers to A^f; — > oo. This is already obvious 
from the parameterisation HOI yielding Q = I on all scales for 
vy = - 1 . In this case, the sound speed becomes irrelevant and dark 
energy clustering cannot be detected. 



3 3D WEAK COSMIC SHEAR 
3.1 Convergence 

In the presence of a gravitational lens, the observed image points 
6 of a galaxy differ from their true positions j8. In a locally linear 
approximation, the mapping S i-» j8 is described by a matrix 



l-K-Ji 
-72 



1 



-72 

- K + yi 



(14) 



The convergence k determines the magnification of the image, the 
shear (71,72) determines its deformation. Both convergence and 
shear are calculated from second derivatives of the lensing potential 
</>, e.g.. 



1 

2 A. 



(15) 



The lensing potential is a projection of the Newtonian gravita- 
tional potential O. In a flat universe. 



Jo 



XX' 



(16) 



vihexe. X , x' denote comoving coordinates. The convergence k thus 
depends on the gravitational potential along the line of sight, which, 
in turn, is given by the density fluctuations. In this way, gravi- 
tational lensing can be used to probe the density field, without 
relying on galaxy bias models i ljain & Selia5ll997l : IHu & White! 
2001.). For general treatments o f weak grav i tationa l lensing, see 
Bartelmann & Schneidej j200lh : lBartelmamj j20I0bl) . 

In a region of the sky covered by a weak lensing survey, the in- 
dividual convergences k and shears 7; of the galaxies together allow 
to study the two-dimensional fields K(&,(fi), ip). If the galax- 
ies' distances are known (e. g. by a photometric redshift measure- 
ment), the fields become three-dimensional, k(x, &, f), Jiix^ V)- 
3d weak cosmic shear is a means to study the statistical properties 
of these fields (Heavens 2003; Castro et al. 2005). 

The importance of the three-dimensional information for 
weak lensing precision tests of str ucture formation has first been 
studied for tomography jHull2002ah . Also the use of spectroscopic 
redshifts instead of a photometric method has been considered 
(Ishak&Hirata2005). 

The statistics of the convergence field k are hardly directly 
observable. But since the statistics of convergence and shear are 
equivalent, we may use the convergence k instead of 7, in our theo- 
retical calculations. 

The first step in a 3d weak cosmic shear calculation is a com- 
bined Fourier and spherical harmonic transform, x ^ k, (&, tp) 
(1, m). 
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Kgm(k) ■■ 



(17) 



By means of the growth factor g{k, a), we may replace the grav- 
itational potential in Eq. l |16t by the potential of today, (b(„,(k) = 
g(k, a) (t>^l^^^{k)/a. Note that in clustering dark energy scenarios, the 
growth factor is scale-dependent. In the transformed variables, 
Eqs. l |15b , l |16l ), and ([7} take the simple forms 



l({ + 1) 



Uk,k')<S'l^{k'), 



(18) 

(19) 
(20) 



where we have, foUowing lHeavenj i2003h . introduced the quantity 

Ut{k,k')=- f /dxMkx) r dx'^^^Mk'x')^^^^i2l) 
n Jo Jo XX « 



and used the summation convention 

A{k,k')B{k',k")= f k'-dk' A{k,k') B{k' ,k"). 
Jo 



(22) 



The appearance of Q on the right-hand side of the Poisson 
equation accounts for the direct contribution, oc p^^A*"^', of dark 
energy perturbations to the gravitational potential. The indirect and 
dynamical effect of dark energy clustering on the evolution of mat- 
ter perturbations p„,A'"'' is accounted for by the modified growth 
index y as outlined in Sec. 12. 31 



3.2 Estimator 



In Sec. 13.11 we have seen that the convergence Kc„(k) is inti- 
mately connected to the density fluctuation field ^(,^{k) by virtue 
of Eqs. ( list to MQ\ . In other words, we can, e. g., use the conver- 
genc e {Kc„,(k) K{'„,'(k ')) to probe the matter power spectrum P^ik). 

iHeavensI (12003) has shown how to construct an appropriate es- 
timator for a weak lensing survey including photometric redshifts. 
The two main ingredients of this estimator are: 

(i) The inclusion of the uncertainty of the redshift measurement, 
for simplicity assumed to be a Gaussian with width cr- equal for all 
galaxies. 



p(x'k)dx' 



1 



^/2no 



■ exp 



(z-z'y 



20-2 



dz'. 



(23) 



We use a typical figure of cr. = 0.02 (Heavens"2003'). An extension 
of the formalism allo wing for individual redshift errors is possible 
( lKitchingetai]|201ll) . 

(ii) The survey's galaxy distribution encoded in the number den- 
sity n(x) = n(x) assumed rotationally symmetric. It constitutes a 
statistical weight favoring distances x where the density of galax- 
ies is higher. We use the forecasted shape 



n(z) dz oc z~ exp 







- 


m 







dz 



(24) 



for the Euclid survey. Here, we assume 100 galaxies per square ar- 
cminute, zo = 0.64, and = 3/2, yi elding a median redshift of 
= 0.9 (|a mara & Refregiedl2007h . For convenience, we con- 
sider the idealised case that the full sky is covered. For a realistic 



sky coverage fs^y < 1, the errors scale approximately by /, 



-1/2 
sky 



We may then define the estimator k[„, for the convergence K[„, in 
terms of the actual convergences Kg of galaxies g as the harmonic 
transform 



Rem(k) ■■ 



(25) 



galaxies g 



As explained above, the cosmic shear would be better suited for the 
analysis of actual observational data. The expectation value of kf^ 
is 

ke^ik) = Z[{k, k') Mtik', k") KUk"), (26) 
with the summation convention Ml\ and the quantities 

Zi{k,k') = ^jx"dx'jdxp(x'\x)Mk'x)jf(kx'\ (27) 

Me(k,k') = Ij X-dxjf{kx)jiik'x)n{x), (28) 

taking account for the two main ingredients stated above. Abbre- 
viating the product Bt{k,k') = Zi(k,k") Miik" ,k"')r}[{k"' ,k'), the 
covariance of k[„, in terms of the matter power spectrum P^{k) reads 

S({k,k') = {ki,n{k)kiJk')) 

= A' B,(kX)^-^^^^^^^P^Bdk',k'') (29) 

withA = ^4;rGpO. 

The full covariance Cf(k, k') is obtained by adding the shot 

noise N[(k,k') = ^Mi(k,k') with crl = 0.1 l lHeaven j l2003h . 
This neglects the non-zero correlation between the ellipticities of 
neighbouring galaxies due to intrinsic alignments ( Heave ns et al.l 
[2OOQ ; S chafe r 2009). This small-scale elfect, however, does not 
afi^ect our analysis of the large-scale consequences of a dark en- 
ergy speed of soun d. Further systematic efi"ects have been studied 
dMarch etaLll201ll) but, in gen eral, do not seem to have a strong 
impa ct on parame ter estimation dHuterer et alj200^ ; lKitching et al.l 
i200a lTakad"a^Jaini2009i) . 



4 PARAMETER ESTIMATION 

We apply a standard Fisher information matrix method to inves- 
tigate possible future parameter constraints from upcoming weak 
lensing surveys. The Fisher information matrix F^y is a square 
matrix whose indices label (cosmological) parameters Pij,py. We 
choose the parameters p^,Py 6 j£l„, A^, /i, n,, w, logj^ cj), assum- 
ing flatness: Q.de = 1 - fim- The Fisher matrix determines strin- 
gent bounds on how precise a parameter p^ can be constrained. If 
all parameters are estimated from the experimental data, the indi- 
vidual uncertainty Ap^ does not go below th e Cramer-Rao bound , 
^Pfi ^ -\/(F"')«i (for an introduction, see iTegmark et al.lll997h . 
The Cramer-Rao bound not only applies to individual parameters, 
but it also determines optimal confidence regions for a set of param- 
eters. For two parameters p^ and py, the corresponding coefficients 
of F"' are a quadratic form defining an error ellipse. 

Formally, the Fisher matrix is defined via the likelihood L, 



F,y = {-d,dy\nL) 



(30) 



The likelihood L = L(k(„\p) is the probability for an experiment to 
measure the value ke,,, for the estimator given cosmological param- 
eters p. 

The cosmological parameters enter the likelihood in two ways. 
First, they predict a power spectrum f^ik) and a growth function 
g(k, a), which, by Eq. i29\ , are decisive quantities for the covari- 
ance of the estimator. Second, they define the background evolution 
and hence the distance measures entering the quantities Zf, M/, and 
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If the likelihood L is a multivariate Gaussian in the data with 
covariance matrix C, the Fisher matrix is given by 



F^, = itr[c-'(a^C)C-'(5,C)], 



(31) 



sensitive to the derivatives of the covariance C with respect to the 
cosmological parameters. 

For our estimator k[„,(k), the covariance matrix carries the in- 
dices (f,m,k,k'). Since different modes I and m are uncorrelated, 
the covariance matrix C splits into blocks. Further, the covariance 
Ci{k, k') from Sec. |3.2| is assumed to be independent of m whereby 
all 21 + \ blocks for a given { are identical. In terms of the covari- 
ance C({k, k'), we may reformulate Eq. Olt to 



(32) 



strument properties as listed in Table 1 of lHollenstein et al.l ( l2009h . 
namely a sky coverage f^iy = 0.65, a beam width Sfwhm = 7 ar- 
cmin, temperature noise = 28 yuK arcmin, and polarisation noise 
Ae = 57 fiK arcmin. For the numerical calculation of the theoreti- 
cally predicted multipoles, we employ CAMB. 



5 METHOD 

In principle, we have already collected the ingredients for our 3d 
weak lensing calculations, namely the covariance of the estimator, 
cf. Eq. J29b , and the Fisher information matrix, Eq. i32i . Due to the 
presence of multiple nested integrals, the actual calculation is in- 
volved and motivates the choice of adequate numerical approaches 
and techniques. We present our strategies in this section. 



It should be kept in mind that the Cramer-Rao bounds are re- 
alistic estimates of the actual constraints only if the likelihood L 
as a function of the parameters p is a Gaussian. This is often vi- 
olated in the case of parameters that are difficult to measure and 
therefore weakly constrained, such as the sound speed parameter 
Cj. The broad likelihood extends to regions where the dependence 
of the matter power spectrum P{k) on c^, cf. Fig. [T] cannot be ap- 
proximated linearly ( iBallesteros & Lesgourg ues.,2010.). This also 
affects the weak lensing convergence spectrum considered in this 
work. Figure [T] suggests that the logarithm logjgc^ is a more nat- 
ural parameter to describe the reaction of the model to variations 
in the dark energy speed of sound. We thus choose log[Q as a 
model parameter in our analysis but emphasize that the Cramer- 
Rao bounds we calculate are only rough estimates of the actual 
future constraints. 

A very practical feature of the Fisher matrix is its additivity. 
Given Fisher matrices F^"^', F^^' for two independent experiments 
A and B, the joint Fisher matrix providing the combined parameter 
constraints is simply F^^^^' = F),'^' + F^^'. This follows directly 
from the multiplication of the corresponding likelihoods and the 
definition of the Fisher matrix, Eq. i30\ . In our case, we can use 
this formalism to include prior information from other experiments 
than weak gravitational lensing. 

As prior information, we use a Fisher matrix F^™^' for the 
cosmic microwave background (CMB) based on forecasts for the 
Planck satellite. We include temperature (TT), polarisation (EE), 
and the cross-correlatio n spectrum (T E). W e calculate the Fisher 
matrix 

^(CMB) following IPerotto et alj jioO^. The predicted noisy 
spectra C^'^ (temperature only), C^^ (F-mode polarisation), and 
Cj^^ (cross-correlation) are encoded in a 3 x 3 matrix. 



(2^+l)/Ay 



1^1 



■TE qTT 



f'TE f^TT 
f^TE f^EE 



TE /^EE 



(33) 



with a fraction /jj^y of the CMB covered. From this, we evaluate the 
Fisher matrix. 



5.1 The quantities Z(, M[, and ri( 

The expectation value k[„, of the 3d convergence estimator, 
Eq. l l26t . mainly is the application of Z[ l l27t . M( i28t . and rjc Ml\ 
on today's gravitational potential <I)^^(fc), 

ki,„ oc Zi{k,k')Mi(k\k")rji{k",k"')<S'l,(k"'X (35) 

where each multiplication corresponds to a integration accord- 
ing to the convention, Eq. l l22t . We have introduced the shorthand 
Bi(k, k'") = Z[(k, k')M[{k',k")T]({k", k'") for the product. 

Before we explain an elaborate way to calculate Bf with high 
precision, we first turn to a simplified approximate approach. Re- 
calling that the sequence of functions 



Mx)= v V|jv(a ^ = ><x+l), 

approaches the Dirac delta function 6o for t 
sufficiently large t, use the approximation 



1 



jrikx) ■' 



(36) 
oo, we may, for 

(37) 



In this approximation, the quantities Zf, Mf, and r]f take simple 
forms, namely 



Mc(k,k') « Lni^l^6D(k-k'), 
rit(k,k') « 



k' -ks{k',a(y/k')) 



for k ^ k' else. 



(38) 
(39) 
(40) 



k^k' a(y/k') 

Calculating the final product B( now does no longer pose difficul- 
ties. 



B[(k,k')^ 



_y_ 

k^k 



2g(k',aiy/k')) n ly \ ..X-jr 
- aiylk') I'^Ak^)"^^!^- ^''^ 



We compare this approximate result with the full expression in 

Fig.m 

Although useful for a first impression, these approximate re- 
sults do not allow for a precision calculation of the covariance 
Cc{k, k'). We thus develop a more sophisticated strategy. 



f=2 PP'.QQ' 

with the indices PP', QQ e {TT, EE, TE]. 

Our forecast bases on expected pro perties of the Planck satel- 
Ute jHoUenstein et alj200^ : lKnoxll99"3) . We adopt the expected in- 



5.2 Covariance 

The signal and noise parts of the covariance matrix are given in 
Sec. 13.21 While the noise part oc M{ is uncomplicated, the di- 
rect evaluation of the signal S c J29b would, in a first step, require 
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wave vector k [Mpc]' 



wave vector k [Mpc]' 



Figure 2. The matrix B((k,k') = Zr{k,k") Me{k" ,k"')ric{k"' ,k') for e = 
10 (upper surface) with the difference between the full integration and the 
approximation given by Eq. j4U (lower surface and contours), which shows 
a small oscillatory feature close to the steep edge of Be(k, k') amounting to 
less than 10% of the amplitude. 



the calculation of Z(, M(, and rjt, which contain highly oscillat- 
ing integrands, cf. Eqs. J28l >. Olb . In a second step, the prod- 
uct BiiKk') = Zc{k,k")M[(k",k"')T][(k"',k') has to be calculated. 
Taken together, these are seven nested integrals. Calculating the 
signal covariance S [ then requires two further integrations. 

Fortunately, the orthogonality relation for spherical Bessel 
functions. 



k^dkj,(kx)jc(kx') = ^6o(x-x'X 



(42) 



r 

can be used to solve several k integrals analytically. The remaining 
expression for Bg reads 

Be(k,k') = 1 J x'dx'jdkx') j P(x'\x)n(x) feik' ,x), (43) 
where 

J"^" , , X - x' s(k, a') 
dx'Mkx') ^ , ' . (44) 

XX a' 

The number of nested integrals in the calculation of B[ is reduced 
to three. 

We will show that an efficient evaluation of the inner integral 
in Eq. i43t is possible using a Fast Fourier Transform (FFT). We 
therefore have to sample f( at discrete coordinates [Xj]- For each 
we need not calculate the full integral j44b but only an integral 
from^j_i loxj - This is possible once we write the integral in a way 
that makes the integrand independent of the integral bound, 

fc(k,x)= —Jeikx) dx je(kx) ; — -(45) 

Jo X a' X Jo 



a' 



In redshift space, the conditional probability p(z'\z) is a Gaussian, 
cf. Eq. J23l >. Inserting this property, we reformulate the inner inte- 
gral in Eq. i43i as a convolution. 



r 



'kp(x'\x)"{x)Mk',x) 



d£ r 

dx' Jo 



dx. 



dzpiz -z)x\^n(x(z))Mk',z) 



(46) 



For convolution integrals, fast solving methods exist. This is due to 
the convolution theorem stating that the Fourier coefficients of the 
individual functions can be multiplied to give the Fourier coeffi- 
cients of the convolution. The Fourier transform of the Gaussian is 
again a Gaussian and thus analytically known. In the last factor, we 
use the sampled values of fc to perform a Fast Fourier Transform. 



5.3 Fisher matrix 

Once the covariance Ce(k,k') is known, the Fisher matrix F^,, can, 
in principle, be calculated according to Eq. J32| l. In terms of linear 
algebra, the Fisher matrix is given by a trace, which is a basis- 
independent operation. This opens the possibility of calculating the 
covariance in another basis, allowing for a more efficient numerical 
calculation. 

Working with the tools of linear algebra, we find it more trans- 
parent to abandon the summation convention I l22t for a moment and 
to work with standard notation instead. All earlier expressions can 
easily be reproduced if quantities of the type A[(k, kf) are replaced 
by ordinary matrices 



A^j., = ■\fk^Ac(k,k') VFm 



(47) 



with a discrete step size Ak. The additional factors automatically 
reproduce the summation convention once a matrix multiplication 
is performed, 2^. k^Ak — > J k^dk. 

Let us search for an orthogonal transformation of the co- 
variance matrix C^j,, = VFm Ci(k,k') VP^aI, 



c' = (t^)' C'Y. 



(48) 



A good choice would, when applied on C'^, produce the orthogonal- 
ity relation for spherical Bessel functions. Such a choice is given by 



^^BAkjtikp) ^Ip^p. 
The noise part N'^ oc becomes particularly simple. 



= Z'^*p'^»''^*'p' = "(p)v- 

k,k' 



(49) 



(50) 



For the transformed signal part, s'^ = (t'^ T'^, the product B'^ is 
transformed from the left-hand side only. 



^Pk' ~ XI ■ 



(51) 



In fact, this transformation further simplifies by virtue of the 
orthogonality relation. 



< I dx' jeik'x' 
Jo 



dxp(pk)n(x) X 



X-X' g(k',a') 
XX' a' 



ylk'^Ak. 



(52) 



Applying the matrix T^^ introduced above on a quantity A^^., 
can be understood as undoing the transformation x ^ k in the 
harmonic transform, cf. Eq. (17}. This means that the Fourier mode 
k is replaced by a comoving distance, now labelled by p. Hence, the 
application of T^^ avoids unnecessary integrations originating from 
the harmonic transform. 

Finally, we have all the necessary tools for an efficient calcu- 
lation of the Fisher matrix at our disposal. 
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6 RESULTS 



The Fisher matrix formalism, cf. Sec. [4] and our numerical meth- 
ods, explained in Sec. [5] enable us to estimate which constraints on 
the dark energy sound speed will be possible with the weak lensing 
data of Euclid. The constraints depend, however, on the assumed 
fiducial parameters since the Fisher matrix is defined by derivatives 
at these points l |31| l. Unfortunately, the dependence of sound speed 
constraints on the fiducial values of both the sound speed Cj itself 
and the equation of state w is very strong. This is illustrated by the 
scale /ieff introduced in Sec. 12. 31 Eq. ( 1131 ), below which dark energy 
clustering is not expected to be observable. This scale is a function 
of both c~ and w, it exceeds the Hubble horizon for » 1 + w. In 
particular, the most natural fiducial value for w mimicking the stan- 
dard ACDM model, i.e., w — » -1, is a singular choice, A^g — » oo. 
The question of how well the sound speed can be constrained 
crucially depends on how close the equation of state w is to the 
value - 1 . 

In order to explore this behaviour quantitatively, we apply 
the Fisher matrix formalism to estimate the uncertainties of the 
dark energy sound speed and equation of state as functions of the 
fiducial values and w. In Sec. O we argued that a natural pa- 
rameter to constrain is the order of magnitude logm rather than 
itself. The relative error on the sound speed approximately is 
ACj/Cj a; In(IO) Alogm c~. This becomes imprecise for large un- 
certainties. In the case of the equation of state, we estimate Aw/|w|. 
For simplicity, we assume all other cosmological parameters to be 
exact ly known, fixed to the WMAP7 recommended ACDM param- 



tiy 

eters dKomatsu et alj20Tlh . 

The Fisher matrix is then a 2 x 2 matrix, and the uncer- 
tainties are estimated as explained in Sec. [4] We combine CMB 
and 3d weak lensing constraints. In our numerical calculation, the 
multipoles £ run from /'min = 2 to /'^ax = 50, the mode k from 
kmm = iO'^ Mpc"' to ytraax = 10"' Mpc"' in A'* = 200 equidistant 



steps. The included redshift range is Zm 



W-" to Zn, 



10 in 



A'; = 1000 steps. For the CMB Fisher matrix, we include, as in all 
subsequent calculations, multipoles from ^ = 2 to ^ = 2250. The 
uncertainties are shown in Fig. [3] These results should be taken as 
a first approximation due to the limitations of the Fisher formalism 
when applied to weakly constrained parameters, cf. Sec.|4] 

The lower figure, which shows the relative error Aw/|w| on 
the dark energy equation of state w, is easily interpreted. The con- 
straints on the equation of state parameter w are largely independent 
of the assumed sound speed c^. So, at least, the uncertainty in the 
sound speed does not worsen the accuracy with which w can be 
known, nor will a wrong assumption on introduce a significant 
bias on the estimate of w. 

In the upper figure, we see that, conversely, the sound speed 
constraints heavily depend on the fiducial values, as explained 
above. For w > -0.95 and sufficiently small the estimated error 
AlogiQC^ is smaller than one. We may thus hope that the combi- 
nation of 3d weak cosmic shear and the CMB will determine the 
order of magnitude of c^. 

For subsequent calculations, we choose the exemplary fidu- 
cial value = 10"~. For w, the most natural choice, w = -1, is 
not adequate. If we still chose w close to -1, e.g., w = -0.99 or 
w = -0.9, all results would strongly depend on the exact value cho- 
sen. Instead, we decide to go further away from the observationally 
preferred value and use w = -0.8 for illustration. 

Adopting these choices for and w as the fiducial val- 
ues, together with the ACDM WMAP7 recommended parameters 
dKomatsu et alj|20lll) . we now calculate the full Fisher matrices 




-0.95 -0.9 -0.85 -O.S -0.75 -0.7 -0.65 -0.6 
dark energy equation of state w 



-0,85 -0,8 -0,75 -0.7 -0.65 -0.6 
dark energy equation oi slate w 



Figure 3. Cramer-Rao bounds on AlogigCj (upper figure) and Aw/|w| 
(lower figure) for varying fiducial values w and log|Q . The three red lines 
mark (from top to bottom) the scales A^g = 10"", n = 1,2,3 times the 
Hubble radius "H"' according to Eq. )13t in the w-c, plane. 



for our six cosmological parameters: fractional matter density n,„, 
scalar initial perturbation amplitude A,, Hubble parameter h, scalar 
spectral index iis, equation of state w, and sound speed logj^Cj. 
We choose higher numerical precision, t^^^ = 300 and Nj, = 500, 
and avoid nonlinear scales. The resulting confidence regions for 3d 
weak cosmic shear alone and for the combined constraints with the 
CMB are seen in Fig.|4] 

3d weak cosmic shear obviously provides interesting con- 
straints on all the six cosmological parameters included in our anal- 
ysis. Some constraints considerably improve when the CMB Fisher 
matrix is added. This is not true for the sound speed. We empha- 
size, of course, that our choice of fiducial parameters of w and 
is only illustrative. In more realistic cases w a; -1, the constraints 
will be much weaker, cf. Fig. [3] 

The Fisher matrix F^y for 3d weak lensing is obtained from 
a summation of all multipoles £, cf. Eq. i32\ . It is instructive to 
examine which multipoles most contribute to the parameter con- 
straints. We therefore plot the uncertainties of all the parameters as 
functions of the maximum multipole ^niax in Eq. l |32| l, see Fig. [5] 

Let us first consider the parameters other than c^. These show 
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0.28 
; 0.275 
0.27 
0.265 



0.69 0.7 0.71 




-0.84 -0.8 -0.76 
" ^.78 
g -0. 





-2 

iog,„(c) 



0.69 0.7 071 



-0.84 -0.8 -0.76 
" ^.78 
g -0. 



-2 
l<>9,„«:) 



2.4 

2.2 
2 
1.8 
1.6 
1.4 
1.2 
1 

0.8 
0.6 



10-" 




10 

maximum multipole 



100 



Figure 5. Marginalised uncertainties AlogjQCj on the dark energy sound 
speed (top panel) and hp^lpi^ of the other cosmological parameters 
(bottom panel), depending on the maximum multipole /"i^ax from 3d weak 



lensing with a CMB prior Fiducial parameters: w = -0.8, 
CMB prior generally includes all multipoles from 2 to 2250. 



10"^. The 



Figure 4. Estimated confidence elHpses (lo", 2<t, and 3o") for 3d weak lens- 
ing alone (upper figure) and combined with CMB prior (lower figure). The 
constraints on Aj have been rescaled by a factor of lO'. 



two distinct behaviours. The constraints on the parameters w, fl„,, 
and h are strongly improved by going to larger multipoles. The two 
parameters Aj and characterising the primordial scalar pertur- 
bation spectrum are already tightly constrained for low multipoles. 
This is linked to the different sensitivities of the two independent 
observations, 3d weak lensing and the CMB, on these parameters. 
Looking again at the error ellipses in Fig. [4] we see that the con- 
straints of 3d weak lensing alone on w, Q.„„ and h are not much 
weaker than the combined ones. Here, 3d weak lensing can estab- 
lish strong constraints with increasing ^,nax- On the other hand, the 
CMB is more sensitive to A, and n,, whereby 3d weak lensing, 
regardless of cannot contribute very much to the constraints. 

The case of the dark energy sound speed is different. The 
fact that the uncertainty does not decrease significantly with in- 
creasing ^max <: 20 is mainly the consequence of clustering dark 
energy being a large-scale phenomenon, cf. Sec. 12.31 and Fig.[T] 
Plotting the covariances Ci(k,k') of the estimator, see Sec. 13.21 for 
increasing multipoles £, we see how the maximal sensitivity moves 
to smaller scales, Fig.[6l In fact, for low multipoles d, 3d weak shear 
probes the scales of interest where dark energy clustering mainly 
occurs. The maxima seen in Fig. [6] are related to the fact that the 



galaxy distribution n(x), Eq. l l24t , peaks at a comoving distance 
X, characterising the survey. Approximating Bessel functions by 
Dirac deltas, cf. Eq. l l37t , this distance roughly corresponds to the 
scale k X tlXf This explains the shift of the maximum for varying 
E observed in Fig. [6] 

Another way to study the { dependence for the parameter con- 
straints is to look at the direct contribution of a multipole £ to the 
diagonal elements F^^ of the Fisher matrix, Eq. l |32| l. These quan- 
tities can be interpreted as a (squared) sensitivity .s^ per £ mode, 

i;;Htr(C7'<9„Cf)' (53) 

Another way of interpreting Eq. i53t is that describes the deriva- 
tive of the measurement with respect to a cosmological parameter 
normalised by the noise of the measurement, such that it assumes 
large values for strong dependences of the signal on the cosmologi- 
cal model and small noise contributions. At the same time, s'^ is the 
contribution to the Fisher matrix entry for the parameter p,, from 
each mode kc, of the convergence field. The number of modes for 
each £ is given hy 2£ + 1. 

We show the sensitivity s'^ for the cosmological parameters 
in Fig. |7] Again, we observe that the sensitivity of 3d weak cos- 
mic shear on the dark energy sound speed c] mostly comes from 
the first multipoles. This confirms the impression already obtained 
from the covariances C((k, k') shown in Fig.[6] and emphasises the 
fact that the influence of a nontrivial sound speed is a large-scale 
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Figure 6. Qualitative results for the covariance C((k,k') for ^ = 5, 10,20 and w = -0.8. Brighter regions mark larger values. The red dashed lines mark the 
scales k^fi = 1 /A^g in the case cj = 1 . 



10' 



10' 



10" 



10-' 



10"" 



10"'= 



10" 



10 

multipole I' 



100 



Figure 7. Sensitivity .v^ as defined by Eq. )53t for the cosmological param- 
eters. Fiducial values: w = -0.8, c~ = 10"^. 



phenomenon. At the same time, the plot explains the two orders of 
magnitude difference in constraints on the dark energy sound speed 
compared to the other cosmological parameters. 



7 SUMMARY 

We have studied the potential of 3d weak cosmic shear to constrain 
a possible clustering of dark energy with the data of next generation 
surveys. We have parametrised the clustering dark energy compo- 
nent by two parameters characteristic for a generic cosmological 
fluid, its equation of state w and its (rest-frame) sound speed c^. 

(i) For the 3d weak cosmic shear analysis, we have developed 
adequate numerical tools allowing for an efficient calculation of 
the covariance and Fisher matrices. These tools were shown to be 
numerically very efficient, which ultimately allowed us to sweep 
through the dark energy parameter space, while retaining sufficient 
numerical accuracy. 

(ii) The capability of future observations to constrain c^, and 
thereby the clustering of dark energy, strongly depends on the dark 
energy equation of state w. If w is close to -1, dark energy per- 
turbations are mainly present at very large scales possibly outside 
the Hubble horizon. The effects of clustering dark energy would 



hardly be observable at all if » 1 -I- w. As, indeed, current obser - 
vations prefer values of w very close to - 1 jKomatsu et al.ll2oTlh . 
this seems to be the decisive caveat. 

(iii) Due to the sensitivity of sound speed constraints to the 
assumed exact value of w, we have estimated the uncertainties 
A logiQ c J and Aw/|w| as functions of the fiducial parameters (w, cj), 
cf. Fig. [3] The numbers are combined constraints based on as- 
sumed properties of Euclid (JTeaven s 2003) and the Plan ck satellite 
jPerotto et ai]|2006l ; lHoilenstein et aLir2009l ; lKnoxlll995h . 

(iv) For the considered range of fiducial parameters (-0.99 < 
w < -0.6, lO""* :£ < 1), the estimated constraints on the sound 
speed AlogiflCj vary between the extreme cases of 0.1 and 3. If 
w > -0.95, the combination of Euclid and Planck is promising to 
constrain the order of magnitude of c J provided that the true sound 
speed is small enough. This would be considerable progress com- 
pared not_onl^;_to_£onstrants_go current observational 
data dde Putter et ai]|2010l : iLi & Xiall2010l) but also, for most of 
the parameter space, to the constraints e xpected from weak lensing 
tomography and galaxy surveys alone jSapone et al.ll20l3) . Con- 
straining Cj within one or two orders of magnitude could also be 
possi b le with Planck and next-generatio n galaxy surveys (Ta kadal 
|2006| . iBallesteros & LesgourguesI I2OICII) o r for neutral hydrogen 
surveys ( Torres-Rodriguez & CressI 20071 : iTorres-Rodriguez et al.l 
l2008h . Although not our focus here, the constraints of 3d weak 
cosmic shear together with the Planck satellite on the dark en- 
ergy equation of state w are worth mentioning; in fact, according to 
Fig.[3]and Fig. [4] w can be constrained below the p ercent level (for 
other constraints from 3d weak cosmic shear, see lHeavensll2003l: 
iHeavens etai]|2006h . Additionally, the true value of w largely de- 
termines the accuracy on the sound speed Cj. 



There are two very different conclusions one could draw from 
these constraints on c^. First, we may regard = 1 and w a -1 
as the natural values as they refer to unclustered dark energy such 
as a cosmological constant. Then, small deviations from = 1 
are interesting; but these seem hardly observable in next generation 
experiments. Second, however, one may argue that is a com- 
pletely unknown parameter with a natural range from to 1 ; then, 
Euclid and 3d weak lensing could single out an order of magnitude 
in which lies. This could be a decisive step for discriminating 
between different dark energy models. 
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